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The model describing market dynamics after a large financial crash is considered in terms of 
the stochastic differential equation of Ito. Physically, the model presents an overdamped Brownian 
particle moving in the nonstationary one-dimensional potential U under the influence of the variable 
noise intensity, depending on the particle position x. Based on the empirical data the approximate 
estimation of the Kramers-Moyal coefficients Di,2 allow to predicate quite definitely the behavior 
of the potential introduced by D\ — ~dU/dx and the volatility ~ \/Th- ft has been shown that the 
presented model describes well enough the best known empirical facts relative to the large financial 
crash of October 1987. 
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I. INTRODUCTION 

The dynamics of the financial markets has been at- 
tracting attention of the physics community for two 
decades p], y, @, 11] • During this time a large volume 
of the empirical researches has been done. Of specific 
interest is the investigation of the behavior of the market 
in the time periods of the large financial crises when the 
statistical properties of the market are drastically distin- 
guished from its properties in quiet days 3, 5, 6, 7]. The 
empirical analysis has found the occurrence of a number 
of peculiarities in market dynamics appearing in these 
periods. Originally they have been revealed when in- 
vestigating the large financial crash on 19 October 1987 
(Black Monday) at New York Stock Exchange. However 
these peculiarities are not specific for October 1987 and 
have been later revealed in the financial crashes of 1997 
and 1998 years for different countries and Exchanges 7]. 

Firstly, it has been established that the large finan- 
cial crashes are outliers which are not possible within the 
scope of the typical price distributions. For their realiza- 
tion the complementary factors absent in quiet days are 
needed [1,0]. 

Secondary, the emergence of the certain periodic pat- 
terns in the dynamics of the different financial assets are 
detected both before the crash, and immediately after it. 
In particular, a log-periodic oscillations in time evolution 
of the price of an asset appear before the extreme event 
0, H, B E3, HI El- The aftershock period is charac- 
terized by an exponentially decaying sinusoidal behavior 
of a price [7], [8| . The characteristic behavior of volatil- 
ity is also revealed after the crash. In the moment of 
the crash of October 1987 the implied volatility of the 
S&P500 index makes a shock and then sharply decays as 
power law decorated with log-periodic oscillations [3, Q • 



In the day of the maximum drop the central part of the 
empirical return distribution moves toward negative re- 
turns and then begins to oscillate between positive and 
negative returns [6j. 

At last, the relaxation dynamics of an aftercrash pe- 
riod is characterized by what is termed as the Omori law 
which shows that the rate of of the return shocks larger 
than some threshold decays as the power law with expo- 
nent close to 1 for different thresholds pjj, [lj] . 

There are a large number of works devoted to the mod- 
eling of the market dynamics immediately before crash 
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[l6j (see also references in the 
]). These works based on the em- 
pirical observations exploiting the analogy with critical 
phenomena have been directed to the investigation of a 
possibility of predictions of the large financial crash on 
the base of the modeling of the dynamics of different fi- 
nancial indexes (S&P500, DJ, etc). 

The "microscopic" approach has been suggested in the 
work [171 ] where the model of the market dynamics taking 
into account a possibility of occurrence of the crash is 
developed. In physical terms, the model is reduced to 
a motion of a Brownian particle in the stationary cubic 
potential and crashes are considered as rare activated 
events. 

In Refs. [l8|, [lj| the generalization of the above model 
has been suggested for the case of the stochastic volatil- 
ity which was considered within the scope of the modified 
Heston model. Physically, the model represents an over- 
damped Brownian particle moving in the stationary cu- 
bic potential under the influence of the fluctuating noise 
intensity. 

Since the statistic properties of the volatility in the 
Heston model are uniform in a time 20] , both the above 
models describe the uniform stochastic process. In fi- 
nancial terms, this means that statistic properties of the 
market are the same regardless of where market is ei- 
ther in normal regime or close to a crash. It seems this 
is poorly consistent with the assertions that crashes are 
outliers and the statistic properties of the market are 



drastically change in extreme days [a, ufl • While within 
the scope of Brownian motion the assertion that crashes 
are outliers and strong nonstationarity of the market re- 
vealed in extreme days pj] lead to the conclusion that an 
"external" time-dependent potential perturbing the sta- 
tionary potential profile, in which a fictitious Brownian 
particle moves in the case of typical days, is triggered in 
the vicinity of the crash. It is likely the appearance of 
such potential could be traced developing the microscopic 
approach in the spirit of the work [171 ] . It is worth saying 
that the authors of this work have also noted that in the 
vicinity of a crash additional market mechanisms , which 
have not been considered in [17j, must be taken into ac- 
count to remove divergences emerging in their model of 
the aftershock dynamics. 

Another way to detect the influence of the "external" 
nonstationary perturbation is to analyze the empirical 
data concerning the extreme events. It is such approach 
that will be considered in the given work for the modeling 
of the market dynamics during the large financial crash 
of October 1987. 

The paper is organized as follows. In SecHT] in terms 
of the stochastic differential equation (SDE) Ito we con- 
sider the model describing the market dynamics during 
a crash. Sec JIIII is devoted to the analysis of the empiri- 
cal data and the estimation of the Kramers-Moyal coeffi- 
cients. In Sec lIVI the discussion of results and conclusion 
are given. 



II. THE MODEL 
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FIG. 1: The time evolution DJ and NASDAQ indexes. 



volatility can be considered as a function of the return 
x and employ the one-dimensional model §Q with the 
variable diffusion coefficient Dz(x,t) being a measure of 
the volatility. 



As it is known [2l| the coefficients D\ and D% from 
SDE (JTJ) , defining the the stochastic process x(t), coincide 
with the Kramers-Moyal coefficients and are given by the 
equations 

D k (x,t) = lim-M k (x,t,T) (A = 1,2) (3) 

M k (x,t,r)= dx'{x' -x) k P{x' 1 t + T\x,t), (4) 



We consider the model in which immediately before a 
crash and in the short aftercrash period the dynamics of 
a financial assets is described by the SDE Ito: 



dx = D 1 (x,t)dt+ s/D 2 (x,t) dW{t), 
dU(x,t) 



Di{x,t) 



dx 



(i) 

(2) 



where xlt) = lnS(t + l)/S(t) is the daily log-return, S(t) 
is the price of the asset at time t, W(t) is a standard 
Wiener process, Z?i and D 2 are the drift coefficient and 
the diffusion coefficient, respectively. Physically, Eq. ([T]) 
describes a motion of an overdamped Brownian particle 
in a potential U(x,t) [21(, introduced by Eq. ((2]), with a 
variable noise intensity depending on a particle position. 
Like the specific models considered in [13, Us Eq- © 
presents a Markovian process. 

As it has been noted above in the moment of the crash 
of October 1987 the implied volatility of the S&P500 in- 
dex makes a shock with subsequent relaxation to its be- 
foreshock level during several weeks. Simultaneously af- 
ter the index decline, its "stabilization" is approximately 
reached within the same time [3, [8|]. All this allows to 
consider that a peculiar synchronization of their move- 
ments occurs during the crash. In turn, this circumstance 
makes it possible to assume that in this period the asset 



where P(x',t + r\x,t) is the conditional probability dis- 
tribution function of x(t). 

Originally, within a financial context, the method 
of the estimation of the Kramers-Moyal coefficients di- 
rectly from the empirical data has been given in Refs. 
|22l. |23||. Further this method has been applied to differ- 
ent financial assets and exchange rates [2J, [25|, [2s [23 [2S 
[29l . [30l | . In all cases the underlying stochastic process is 
supposed to be stationary (uniform) and the reconstruc- 
tion the drift coefficient D\ from the empirical data gave 
the linear dependence on corresponding return. Accord- 
ing to Eq. ([2]) this leads to the parabolic potential for 
a fictitious particle that is typical for the quiet market 
periods. Here we shall consider the estimation of the co- 
efficients D k for the case of the extreme events of October 
1987 when the market was in a substantially nonstation- 
ary phase p . 



III. THE ESTIMATION OF KRAMERS-MOYAL 
COEFFICIENTS 

"From the opening on October 14, 1987 through the 
market close on October 19, major indexes of market 
valuation in the United States declined by 30 percent or 
more. Furthermore, all major world markets declined 
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FIG. 2: The reconstruction of the potential U. In the figures (a) and (b) the potential surface is viewed from two mutually 
opposite sides 



substantially in the month, which is itself an exceptional 
fact that contrasts with the usual modest correlations of 
returns across countries and the fact that stock markets 
around the world are amazingly diverse in their organiza- 
tion " [3, HI]. In Fig. Q] the time evolution of two financial 
indexes DJ and NASDAQ characterizing different sec- 
tors of USA economics are given. Practically identical 
dynamics of these indexes during the crash is seen from 
Fig. [1] Such behavior is not accidental and is well con- 
sistent with the assertion that large financial crashes are 
caused by appearance of substantial cooperation in the 
behavior of different markets over the world and local 
self- reinforcing imitation between traders [3]. These two 
circumstances led to the substantial synchronization in 
the dynamics of different financial assets simultaneously 
traded both in the common market and in different ones, 
that Fig.l shows. These circumstances allow to suggest 
also that log-returns of different stocks simultaneously 
traded in markets during the crash have the identical 
statistical properties and to consider the probability dis- 
tributions of the whole ensemble of such stocks. Actually 
this approach has been used in Ref. |g|. 



In accordance with the above said we consider the en- 
semble which consists of n = 467 stocks traded in Stock 
Exchanges NYSE, NASDAQ and AMEX in the period 
from 13 October to 13 November 1987 I31J each charac- 
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FIG. 3: The reconstruction of the diffusion coefficient D2. 



terized by the daily log-returns 
Si(t + 1) 



Xi(t) = In- 



Si(t) 



1,2. 



,.n, 



where Si (i) is the close price of i - th asset on day t and 
the time is measured fronr the moment of the crash of 19 







FIG. 4: The potential surface U given by Eq. jTJl 



October. 

Apart from the assumption that during the crash the 
statistic properties all Xi(t) are identical we shall intro- 
duce a fictitious index S*(t) and consider the empirical 
data Xi(t) (i = I, 2, ...n) as the set of "experimental" val- 
ues of the log-returns of this index detected at moment t. 
Further the estimation of the Kramers-Moyal coefficients 
will be performed just for the index S*{t) with log- return 
x{t). 

According to Eq. ^ the calculation of the coefficients 
Dk includes the operation of taking the limit r — > 0. 
For the case of the high-frequency sampling a sufficiently 
exact estimation can be obtained merely as ratio 



On the long time scale with the minimal step r = 1 
day, Eq. ([5]) can be broken down. Nonetheless, the rough 
but, as it will be seen later, giving definite information 
estimate of the coefficients Dk can be also obtained in 
this case at r = I . Actually the approximation is used, 
meaning an extension of the linear dependence on the 
long time scale: 



M k (x,t; 



rD h (x,t) 



(0 < r < 1) (6) 



D k (x,t) 



M k (x,t-,T) 



(r-0) 



(5) 



Thus the estimate of the Kramers-Moyal coefficients 
has been taken from equality Du{x,t) = Mk(x,t,;r = 
I). For the calculation of the conditional moments 
Mk the empirical conditional densities P^^l^i^i) = 
P(x2t2,Xiti)/P(xiti) (i 2 > h) have been first found 
where P{xit2,xit{) and P(xt) are the joint two-point 
probability density and the one-point probability den- 
sity, respectively. Then the numerical integration has 
been performed in Eq. (TJ|. Once the coefficient D\ has 
been found the potential U is determined from Eq. I|2|). 
The results of the calculations of U and L>2 are shown in 
Figs. [Hand E3 

The periodic structure of the potential surface is 
clearly distinguishable in Fig. [5] fn the regions x > 
and x < alternation of wells and hills is observed, the 
sizes of which rapidly decrease with time. The periodic 
structures of these two regions are displaced in relation 
to each other so that at the displacement along x axes a 
well in the range x > transfers to a hill in the range 
x < and conversely. At t = there is the well of the 
largest depth in region x < corresponding to the mo- 
ment of the crash. The behavior of the potential surface 
with indicated features reconstructed from the empirical 
data is well enough reproduced by the smooth surface 
presented in Fig. (U|) and given by the equation: 



U(x,t) =<p(x,t) 



A sine 



-OLX+j3it 



B sin lux sm^LUit + b) e 
I 



t < 

-ax-m t > 



(7) 



where 



,t) = l 



ax^e- 



-T 1*1 



(8) 



The potential U presented by Eq. ([7|) exponentially de- 
cays with time and included trigonometrical functions re- 
produce the observed periodic structure. Function <f>{x, t) 
is largely used for the smooth joint of two surfaces from 
Eq. (JJJ) at t = and rapidly goes to one with the increas- 
ing \t\. The parameters Eqs. ([7]) and ((5} are determined 
by comparison with the reconstructed surface on the base 



of the least-squares method and presented in Table |TJ 



The reconstruction of the diffusion coefficient D2(x, t) 
is shown in Fig. ([3]). Qualitatively the behavior of the 
diffusion coefficient is quite definite, in spite of the occur- 
ring irregularity. At the moment of the crash D% makes a 
shock with subsequent relaxation to its beforecrisis level 
within a month. In Fig. ([5]) such behavior is reproduced 





TABLE I: Parameters of Eq.7 
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TABLE II: Parameters of Eq. © 
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by the smooth surface with the equation 



(9) 



As before the parameters of Eq. ((9J have been defined by 
the least-squares method and given in Table HU 



IV. DISCUSSION AND CONCLUSION 

In previous section we have estimated the Kramers- 
Moyal coefficients D\ and D 2 from the empirical data. In 
spite of the sufficiently rough estimate, the reconstructed 
surfaces for the potential U(x,t) and the diffusion co- 
efficient D 2 , presented in Figs. [2] and El can be quite 
definitely approximated by the smooth surfaces with the 
needed properties. 

Firstly, it has been noted in Introduction that in the 
aftercrash period the center of the empirical return dis- 
tribution oscillates between positive and negative returns 
Q . Evidently it is the periodic structure of the potential 
surface presented in Fig. |4] that provides such behavior. 

Secondly, on the base of Eqs. (JTJ) and |(3J) with U and 
£>2 from Eqs. ([7]) and © the time series of the index 
S*(t) has been generated. Its time evolution is presented 
Fig.[S]by the solid line. It has been obtained by averaging 
over 150 the simulated trajectories yielding the initial 
index decline greater than 25% . As it has been noted 
in the early observations [8[ the aftercrash time behavior 
of the S&p500 index is characterized by an expotentially 
decaing sinusoidal function. In our case the generated 
data fit well enough with the function 

S{t) = (Ai e- ait + A 2 e- a2t ) sm(ut + 7) + ^0 , (10) 

incorporating the superposition of two the exponential 
functions with the different exponent ot\ and a 2 . (the 
dashed line in Fig. [5]). The parameters of Eq. (|10p have 
been defined by the least-squares method and given in 
Table HTT1 

Thirdly, as follows from Eq. © at the moment of the 
crash volatility ~ yHj^ makes a shock and further decays 
to the aftercrisis level as power law with the exponent ap- 
proximately equal to 0.3. The similar behavior has been 
also detected for the implied volatility of the S&P500 
index (with exponent approximately equal to 1.) [8J. 




FIG. 5: The behavior of the diffusion coefficient Di given by 
Eq. ©. 
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FIG. 6: The time evolution of the fictitious index S*(t). The 
solid line has been obtained by averaging over 150 the sim- 
ulated trajectories yielding the initial index decline greater 
than 25% . The fit with Eq. ffTTTf) is shown by the dashed line. 



Finally, we have examined the fulfillment of the Omori 
law within the framework of the given model. To this 
end the cumulative number N(t) of the log-returns in in- 
terval [0, t], whose absolute values exceed a given thresh- 
old value £ has been calculated for different t. As it is 
seen from Fig. [7] the generated data are described well 
enough by power law t 1_n , that is characteristic for the 
Omori law, with fi = 0.626 for the threshold £ — a and 
ft = 0.704 for I = 1.5cr where a = 0.0685 is the standard 
deviation of log-return x(t) computed over the entire pe- 
riod of 25 days. The values of fl are well consistent with 
the data found by Weber at al (l4j . 





TABLE III: Parameters of Eq. (JTUJ) 
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FIG. 7: Cumulative number N(t) of log-returns in interval 
[0,t], whose absolute values exceed a given threshold value £ 
for t — a and 1.5<r, where <r = 0.0685 is standard deviation 
of log-return x(t) computed over the entire period of 25 days. 
The generated data are shown by circles, the solid lines are 
power law fits. 



Thus the presented model, it seems, allows to describe 
the best known empirical facts relative to the aftercrash 
dynamics. To some extent this circumstance can serve as 
validation of performed estimation of the Kramers-Moyal 
coefficients. 



In conclusion, the given work originates from the em- 
pirical observation that the large financial crashes are 
outliers and during the crash the statistic properties of 
a market are drastically changed in comparison with its 
typical properties. Based on these observations the phys- 
ical model presents an overdamped Brownian particle 
with the variable noise intensity denned by the diffusion 
coefficient Z?2, moving in the nonstationary potential U. 
The market mechanism that triggers such potential is 
related to an imitation between traders increased up to 
a critical the point and a panic reigning in a market in 
the days of the crash. On the base of the empirical data 
within the scope of the Markowian approximation the 
Kramers-Moyal coefficients have been estimated with the 
subsequent determination of the potential U(x, t). As it 
has been shown the given model reproduces well enough 
the best known empirical observations detected in the 
days of the large financial crash of October 1987. 
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